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Einstein-Podolsky-Rosen steering is a manifestation of quantum correlations exhibited by quantum systems, 
that allows for entanglement certification when one of the subsystems is not characterized. Detecting steerability 
of quantum states is essential to assess their suitability for quantum information protocols with partially trusted 
devices. We provide a hierarchy of sufficient conditions for the steerability of bipartite quantum states of any 
dimension, including continuous variable states. Previously known steering criteria are recovered as special 
cases of our approach. The proposed method allows us to derive optimal steering witnesses for arbitrary families 
of quantum states, and provides a systematic framework to analytically derive non-linear steering criteria. We 
discuss relevant examples and, in particular, provide an optimal steering witness for a lossy single-photon Bell 
state; the witness can be implemented just by linear optics and homodyne detection, and detects steering with a 
higher loss tolerance than any other known method. Our approach is readily applicable to multipartite steering 
detection and to the characterization of joint measurability. 

PACS numbers: 03.65.Ud, 03.67.Mn, 42.50.Dv 


Introduction — Entanglement certification llj] is a major 
topic within quantum information science, as entanglement 
is behind many applications H. In this task, the knowledge 
of the system’s dimension and the assumption that the mea¬ 
suring devices are trusted, i.e. they operate as prescribed, are 
usually taken for granted; yet a breakdown of any of these 
could undermine the entanglement certification yD. Moti¬ 
vated by quantum cryptography, where minimal assumptions 
are desired, it was realized that nonlocal entangled states, 
i.e. states that violate a Bell inequality 10-0, allow for en¬ 
tanglement certification with systems and devices being com¬ 
pletely uncharacterized 101. Fully-device-independent entan¬ 
glement certification, however, requires the observation of 
Bell violations free of the detection loophole, which is ex¬ 
perimentally very demanding |0 . 

Einstein-Podolsky-Rosen steering iffl] is a type of quan¬ 
tum correlation that is intermediate between entanglement and 
nonlocality ifl^flBll . Witnessing steering in a bipartite system 
implies entanglement certification without any assumption on 
one of the parties, i.e. in a one-sided device-independent man¬ 
ner. Being less stringent than nonlocality, steering is more ro¬ 
bust against experimental noise il2lll4l4l7l] . and its loophole- 
free detection is feasible with present-day technology H. 
Steering also enjoys a plethora of applications, for which plain 
entanglement is not enough while the harder-to-get nonlo¬ 
cality is not required. These applications range from one¬ 
sided device-independent quantum key distribution Cl, ad¬ 
vantages in subchannel discrimination El, secure quantum 
teleportation El, and connections to joint measurability of 
generalized measurements | H 25 ]. 


which give an experimenter the freedom to choose the mea¬ 
surements involved, and allow for an improvement of the de¬ 
tection by performing additional measurements until a vio¬ 
lation is observed lH 34-^. These criteria are based on 


Compared to well-studied entanglement and nonlocality, 
relatively little progress has been achieved about steering de¬ 
tection. A handful of criteria exist Eld, which are how¬ 
ever tailored to specific measurement scenarios. Only very 
recently some constructive steering criteria were introduced. 


the useful methods of semidefinite programming Il37l] . and the 
downside in this case is that, so far, they could only be applied 
to discrete variable (DV) systems with not too high dimension, 
due to computational limitations. It is then clear that there ex¬ 
ists still a gap that needs to be filled about steering detection, 
regarding higher dimensional DV systems and general contin¬ 
uous variable (CV) systems. 

In this Letter, we propose a hierarchy of steering crite¬ 
ria that is directly applicable to bipartite quantum systems of 
any dimension, including the case of infinite-dimensional sys¬ 
tems. Our method avoids the dimension problem by utilizing 
moments of observables instead of dealing with conditional 
states, at variance with previous DV proposals. A system¬ 
atic framework is provided for deriving non-linear steering 
inequalities in an analytical manner. To the best of our knowl¬ 
edge, our proposed method is the first instance of a hierar¬ 
chical family of criteria for quantum steering that is valid for 
any dimension, and shares some similarity in spirit and struc¬ 
ture with the hierarchy of moments by Shchukin and Vogel 
1 3811 for CV entanglement detection, and with the Navascues- 
Pironio-Acm hierarchy for the characterization of nonlo¬ 
cal quantum correlations. We show that our approach provides 
optimal moment-based linear steering witnesses for any cho¬ 
sen states and measurements on both parties, including CV 
ones. Furthermore, various previously proposed steering cri¬ 
teria are retrieved as special cases of our unifying approach, 
while new non-linear criteria are derived. Finally, we con¬ 
sider several examples of both DV and CV states, and show 
that our technique allows to beat the current state-of-the-art 
in steering detection of a lossy single-photon entangled state 
with quadrature measurements 13311 . 
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Steering detection — We consider the entanglement certifi¬ 
cation task in which two distant parties, Alice and Bob, each 
holding one half of a quantum state pab of a bipartite system 
(described by a Hilbert space 'Ha ® 'Hb, where 'Ha,'Hb de¬ 
note the Hilbert spaces of Alice and Bob respectively), want to 
verify that they share entanglement. Additionally to this, we 
impose the constraint that Alice’s system is unknown (i.e., un¬ 
known 'Ha), and her measurement devices cannot be trusted. 
This implies that the measurement outcomes Alice announces 
cannot be assumed to originate from a particular observable 
on some quantum state of known dimension. The usual en¬ 
tanglement criteria in this case are inapplicable and we need 
to consider steering criteria to identify any nonseparability be¬ 
tween the untrusted Alice and the trusted Bob 111 ill . 

In this scenario, Alice performs one out of n unknown mea¬ 
surements (often called ‘inputs’) on her half of pab, labelled 
by X = and with probability p{a\x) gets some out¬ 

come a. In principle, Alice’s measurements are arbitrary, but 
one can restrict the analysis to projective measurements with¬ 
out losing generality, because the ancilla needed for a non- 
projective measurement can always be moved to the definition 
of the local state on Alice’s side. Alice announces the corre¬ 
sponding pair {a, x) to Bob, who then tomographically recon¬ 
structs his conditional local (unnormalized) state which 
is of arbitrary, but known, dimension. Bob’s states are de¬ 
fined so that Tr(cr^i^) = p{a\x). For all possible pairs (a, x), 
Bob thus obtains the set called an ‘assemblage ’ 0 . 

From the assemblage alone, they should judge whether entan¬ 
glement was present between their shared systems. We refer 
to this procedure as a steering test. 

More precisely, based on the observed assemblage, they 
must determine whether there exists a separable model, i.e., a 
separable state = 2,1 on'K^^'TYa, and measure¬ 

ments {Ma\x]x for Alice, that reproduces Bob’s assemblage if 
we allowed for arbitrary Hilbert spaces 'H'^ on Alice. If such a 
model does not exist, then the shared state must be entangled. 
A steering test using a separable state Pab, and measurements 
{Ma\x]x associated to each input, necessarily leads to the fol¬ 
lowing form for Bob’s conditional (unnormalized) states, 

= TrA[(Map(8> ]1 b)pab] = '^qxp{a\x,A)p^, 'ia,x, (1) 

A 

where p{a\x,A) - ^AMa\xp'\} and p{a\x) - Tr lo-^ , .]. Assem¬ 
blages of the form ([T]) are called unsteemble ill ill . One can 
also prove that, given any unsteerable assemblage, there al¬ 
ways exist a separable state and projective measurements for 
Alice that reproduce it. Furthermore Alice’s measurements 
can be assumed to be described by mutually commuting ob¬ 
servables S. Intuitively, this follows from the fact that a 
separable model is ‘classical’ on Alice’s side. Therefore, un¬ 
steerability is equivalent to the existence of such a separable 
model. 

Our approach is based upon the fact that Bob’s conditional 
states, cr^i^, on which the steering test is based, are in gen¬ 
eral hard to obtain experimentally when the set of outcomes 


is large, or even continuous, as Bob would need to do tomog¬ 
raphy for every pair {a, x). To circumvent this problem we 
instead consider the more accessible correlations 

® b;) = ^ P{a, b\x. By) = 2 fl^Tr [cr.pB;], (2) 

a,b a 


between the unknown observables Ax - 2 a <^Ma\x (with x - 
!,...,«) measured by Alice, and some known observables By 
on Hb (with y - m) measured by Bob, with outcomes 
(eigenvalues) b. In Eq. (|2]i, g,T > 0 are integer powers, and 
P{a, b\x. By) is the observed joint probability distribution. In 
what follows we will show how to derive tests for steering, 
based solely upon the observed correlations {(Aj ® B')}. 

Moment matrices — The main tool we will use is a moment 
matrix, defined asakx k matrix F with elements 

Pij-iSjSj}, (3) 


where i,j- I,.. .,k, and each operator S , is some (as-yet un¬ 
specified) product of operators for Alice and Bob. As a simple 
example, if Bob’s system is a qubit, one could choose the set 
S - {1 ® ]l,Ai ® A, A 2 ® T, A 3 ® Z) where Bob’s observables 
X, Y, Z denote the three Pauli operators. 

We first remark that such a moment matrix, when con¬ 
structed from physical observables on quantum states, is 
always positive semidefinite, i.e. F > 0. This follows 
immediately, since for any vector v, with elements v,, 
j)vj = jVj)) > 0. The second cru¬ 

cial property is that if the underlying operators satisfy any 
algebraic properties, then the moment matrix inherits addi¬ 
tional structure in the form of linear constraints. For example, 
if two (hermitian) operators commute, [5,, 5 j] = 0, then the 
corresponding elements of the moment matrix are necessarily 
equal, Fij - (S -S j) = (SjS,) - P y,. As a second example, if 


S]Sj - iSk and = 1, then Fij - {SjSj) - - iF\k. 

In the next section we show that these properties allow us to 
construct a steering test based upon moment matrices. 

Detection method based on the moment matrix — Consider 
a steering test defined by a set of observed correlations (| 2 ]i and 
take any set of operators S involving some unknown opera¬ 
tors on Alice’s untrusted side and known operators on Bob’s 
trusted side. Now consider the unknown moment matrix F 
associate to S defined as in Eq. Q. Some of its matrix ele¬ 
ments however are known as they correspond directly to ob¬ 
servable data in the steering scenario: these include moments 
of the form (l2]i, and moments of the form (Aj ® B), with B an 
arbitrary operator in Bob’s trusted operator algebra MM- 
All the other elements are not directly available, since they 
involve products of Alice’s unknown operators ll4^ . and are 
treated as arbitrary (complex, in general)/ree parameters. 

Our main goal is to check whether the observed data could 
be obtained or not by measurements on a separable state. At 
the level of the moment matrix, assuming that the observables 
A, commute imposes some extra linear constraints between 
the elements of F, as discussed above. Additionally, we can 
also impose other constraints on F given the knowledge of 
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Bob’s operators. The idea of our method then relies on search¬ 
ing for values for the free parameters of the constrained F that 
make it positive semidefinite. If no such values are found, then 
the data are incompatible with a model relying on commuting 
observables on Alice’s side, and consequently no separable 
state could give rise to it. 

More formally, let "R denote a particular simultaneous as¬ 
signment of values to all independent free parameters, and let 
denote the moment matrix for commuting measurement 
operators on Alice’s side dependent on such an assignment. 
Then, steering is witnessed from Tk if the latter cannot be 
made positive semidefinite for any possible assignment R of 
the free parameters, i.e.. 


^ 0, '^R {(Aj ® Zip) demonstrates steering. (4) 

As anticipated, Eq. (HJi is the central result of this Letter. 

The proposed method for investigating steerability through 
moments of observables shows many advantages. First, it is 
valid for bipartite quantum systems of any dimension, be it 
discrete, continuous or even hybrid since everywhere Bob’s 
Hilbert space was assumed arbitrary, while Alice was al¬ 
lowed for an arbitrary (discrete or continuous) set of out¬ 
comes. Second, the condition (HJi serves as an infinite hier¬ 
archy of criteria; one may start with a small set of selected 
operators {Si}, that are chosen at will, and can gradually in¬ 
crease this set by adding more moments to improve steer¬ 
ing detection. In particular, the operators (5,) can be chosen 
from the set S of all strings (products) of operators of Al¬ 
ice’s (unknown) observables, and Bob’s (known) observ¬ 
ables By. This infinite set can naturally be partitioned into 
subsets containing all strings of a given length k. For ex¬ 
ample, with only two operators on each side, = {1 ® 1 ), 
= {Ai®11,A2®1, 1®Bi, 11®B2),‘S'^* = {AiA2®11,A2Ai® 
l,Ai ®Bi,Ai ®B 2 , A 2 ®Bi, A 2 ®B 2 ,1 ®B 1 B 2 ,1 ®B 2 B 1 }, etc. 
Third, checking whether there is any assignment of unknown 
parameters which makes a matrix positive semidefinite sub¬ 
ject to linear constraints is an instance of a semidefinite pro¬ 
gram (SDP) which can be efficiently solved for many cases of 
interest. Moreover, the duality theory of SDPs allows us to 
extract linear inequalities which act as witnesses for steering. 

Examples — In the following we consider various families 
of quantum states, and show that the proposed hierarchy gen¬ 
eralizes and includes known steering criteria as special cases. 

(i) 2x2 Werner states — Consider the class of discrete vari¬ 
able two-qubit Werner states H , Pab(w) ^ w\iJ/ )ab{^ I + 
(1 - w)llAij/4, where (|01>ab - |10>ab) is the sin¬ 

glet. To check their steerability, we construct the moment ma¬ 
trix (O defined by the previously mentioned set of observables 
S - {1 ® 1, Ai ® X, A 2 ® T, A 3 ® Z) for Alice and Bob; 




1 <Ai®X) 

<Ai®X) {A\®X^) 
<A2®T) {A2Ai®YX) 
<A3®Z) <A3Ai®ZX) 


<A 2 ® Y) 
{AiA2®XY) 
{A\ ® t2) 
<A 3 A 2 ®ZT) 


(As ® Z) ' 
<AiA 3 ®XZ> 
(AsAs ® TZ) 
<A 2 ®z 2 ) , 
(5) 


Consider the statistics of Alice’s unknown measurements 
Ai,A 2 ,A 3 to originate from spin-measurements X,Y,Z, re¬ 
spectively, on her share of pab- We observe that (Aj ® B) - 
{X'^ ® B}p^g(„',, for k = 1,2 and arbitrary B, and similarly for 
the observable elements that contain A* and A*. Furthermore, 
the commutativity requirement on Alice’s side, together with 
the algebra of operators on Bob’s side (e.g. (A 1 A 2 ® XY) - 
-(A 2 A 1 ® YX)), reduces the number of independent free pa¬ 
rameters to three. One can then numerically check the posi¬ 
tivity of the moment matrix and find that ^ 0, 'iR, for all 
w > Wmin = 1 / V3, which is known to be the threshold value 
for steering when Alice has exactly three inputs lE^ . as is the 
case here. The dual of the SDP gives the following optimal 
steering witness, for this family of states and measurements. 


<Ai ®X) -I- <A2 ® T> -h <A3 ® Z) > - V3, 


(6) 


which is violated by all Werner states with w > 1 / VS, while 
satisfied by all unsteerable states 0 ]. The steering criterion 
(| 6 ]l was derived independently elsewhere ll27ll . and we have 
shown that it is only a special case of our general approach. 

Non-linear criteria can also be derived and, remarkably, in 
an analytical manner. A hermitian matrix is known to be pos¬ 
itive semidefinite ijf all its principal minors are non-negative 
ll45ll . Since T?; is by definition hermitian, > 0 implies 
detF?; > 0 , that can be shown to be satisfied by all unsteer¬ 
able assemblages iff 0 


<Ai ®Xf + {A 2 ® T>^ H- <A3 ® Z)2 < 1. 


(7) 


When applied to Pab(w), steering detection is achieved for w 
down to the known threshold value w^in = 1 / V3- Moreover, 
based on the positivity of the principal minors of Q, other 
non-linear criteria can be derived with two (instead of three) 
dichotomic measurements per site 0 . 

(ii) Gaussian states — Let us first make a small introduc¬ 
tion to Gaussian states |4^ £Z1- A Gaussian state is well- 
known to be fully determined up to local displacements by its 


covariance matrix, defined as ( 74^),7 = Tr [(/?,/?7 + /?y/?,)p^^], 
which contains all second moments for two modes. The vec¬ 
tor R - (qA,PA,qB,PBy conveniently groups the quadrature 
operators qA(B),PA(B) for each mode, satisfying the canoni¬ 
cal commutation relations = i{QAB)jk, where Qab - 

Qa © Gb is the symplectic matrix, with 


£^a — ^b — 


0 1 
-1 0 


To simplify things, using local unitary operations and classical 
communication, we can always bring the covariance matrix 
into the standard form. 


Yab 


A C 
B 


where A = diag(a,a) and B - diag(Z 7 , Z?) are the marginal 
covariance matrices of Alice and Bob and C - diag(ci, C 2 ) 
contains their correlations. 
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We proceed by investigating the steerability of Gaussian 
states in standard form (which implies the steerability of any 
non-Gaussian state with the same second moments thereof) 
using the following set of quadrature observables, = {Ai ® 
ll,A 2 ®l,l®q'B, 1® Pb], while we consider Alice’s unknown 
measurements Ai,A 2 to originate from measurement of the 
quadratures qA, Pa respectively. The corresponding moment 
matrix T Q for Gaussian states in standard form becomes. 




a 

R 

Cl 

0 f 

R 

a 

0 

C2 

Cl 

0 

b 

i 

0 

C2 

—i 

b , 


( 8 ) 


with 'R - (A 1 A 2 ) being the only unobservable free (real) pa¬ 
rameter with commutativity imposed. We can proceed analyt¬ 
ically, by remarking that if were nonsteerable then there 
would exist R such that > 0 which implies detF^ > 0 . 
The latter, is equivalent to dety^s -det A > R^(detB- 1) > 0, 
where for the second inequality we used the property det B > 
1 that all physical states must satisfy 0 . Therefore, all un- 
steerable assemblages necessarily satisfy det jab - det A > 0 , 
while a violation would signal steering since there exist no R 
able to make detF^^j non-negative and consequently F 51 ; posi¬ 
tive semidefinite. The steering condition det jab - detA > 0 
derived here can be shown to be satisfied iff \ 




Jab + *(0a © ^b) ^ 0, 


(9) 


which is precisely Wiseman et al.’s necessary and sufficient 
criterion for the steerabiljW of Gaussian states under Alice’s 
Gaussian measurements 11 111 12611 . Therefore, yet another cri¬ 
terion turns out to be a special case of our approach and this 
time in the CV regime. It is worth remarking that the deriva¬ 
tion of (|9]l presented here made no assumptions about either 
Alice’s uncharacterized system or the Gaussianity of Bob’s 
subsystem (also, see |49]), in contrast to El. 

(iii) Lossy NOON states — Consider now the following class 
of lossy non-Gaussian CV bipartite quantum states, 

pTb =(!-'?) IOO)Aii<00| + 77 |V00V)Aii<A^00V|, (10) 

where \NQ0N)ab - {\N0)ab - \0N)ab) is the well-known 

NOON state useful in quantum metrology if^ . whose imper¬ 
fect preparation is modelled through a mixing with the vac¬ 
uum with probability 77 . For later use, let us define position 
and momentum observables for each party A{B), given N, as 
I^Th _ J_( JN , V N JN) _ '1 

‘7A(B) “ yj2^^A{B) ^ ^A(B)’ \ MB) ^A(B)p 

satisfying [q%lyp‘'^ = i, with [aA(B), 4 (B)] = 1 - 

For N - 1, Eq. ( fTO describes an entangled state produced 
by splitting a single photon (generated with probability 77 ) at 
a 50-50 beam splitter. This state is of theoretical 53] and 
experimental interest | IJ, 3 , and it is very desirable to have 


an experimentally friendly criterion that allows one to certify 
some form of nonlocality in its correlations. To our knowl¬ 
edge, the current best steering detection for p^*^ using only 
quadrature measurements is achieved by a non-linear steering 


inequality proposed by Jones and Wiseman 113311 . which can 
detect steering down to 77 > 0.77 in the limit of Alice having an 
infinite number of inputs, while both Alice and Bob bin their 
outcomes (i.e. for a given outcome a, a value is assigned 0 if 
a < 0, and 1 if a > 0). For comparison, recently proposed en- 
tropic steering criteria OSll . employing (unbinned) quadrature 
measurements for both parties, can be seen to detect steering 
for a weaker 77 > 0.94, while all criteria that involve moments 
of quadrature s up to second order fail to detect any steering 
at all 01 ll |27[ 15511 . We will show that our moment matrix ap¬ 
proach outperforms all the previous methods for these states. 

To make the comparison fair, we also consider that Alice 
only performs two quadrature measurements, but allow Bob 
to measure arbitrary local operators, see Appendix for discus¬ 
sion. To test for steering we use S - {1 ® 1, Aq ® qB,Ao ® 
Pb,A\ ® qB,A\ ® Pb,A^ ® 1 , Aj ® 1,1 ® o’!, 1 ® qsPB, 1 ®> 
psqB, 1 ®> P^}, with the observable data calculated assum¬ 
ing Alice’s unknown measurements A 1 , A 2 are the quadratures 
qA,PA respectively. Here, qA{B),PA{B) correspond to qA(B)’PA(B) 
defined above, with N = 1. The set S defines an 11 x 11 
moment matrix F ©, with two inputs Ai,A 2 associated to 
Alice. Following the steps of the detection method, with the 


observable elements of F computed from the state p^'^ 14211 . 


we employ SDP to efficiently check I©, and manage to detect 
steering for all 77 down to the critical value 77 > | s 77 ^, which 
is lower than what previous methods can achieve. The dual 
of the SDP gives us the optimal linear steering inequality for 
reported in the Appendix, that is violated for all 77 > | 
and satisfied by all unsteerable assemblages. The proposed 
witness involves for Bob local moments of quadratures up to 
fourth order and can be efficiently measured by homodyne de¬ 
tection and linear optics therefore demonstrating the 

experimental feasibility of our proposal. 

For any given N > 1, we can consider the same set S, with 
corresponding observables q^^^gy P^a(B)' tested our 

method up to V = 6 and observed a steering detection down 
to 77 > 77 f'\ with 77 f^ < I (e.g., 77 ® Si 0.61 for N = 6 ). We 
conjecture that steering be detectable with our method for all 
N, although larger values could not be tested due to computa¬ 
tional limitations. We should note however that for A > 1 the 
observables q^^^gy P^a(B) correspond to non-Gaussian measure¬ 
ments that are hard to implement experimentally. On the other 
hand, for A > 1 the feasible quadrature measurements qA{B) 
and pa(B) could not detect steering in the states of Eq. (fTOl) for 
any 77 and for the given set S considered above. 

Conclusion — We proposed an infinite hierarchy of suffi¬ 
cient conditions for bipartite steering applicable to all quan¬ 
tum systems. Other previously known steering criteria were 
shown to be special cases of our approach, both in the discrete 
and continuous variable regimes. An optimal witness for an 
inperfect single-photon entangled state was obtained, which 
was shown to be more resistant to losses than previous pro¬ 
posals, and experimentally accessible with linear optics and 
homodyne detection. In the light of a recently proved equiv¬ 
alence between steering and joint measurability 23.1^1. 
the hierarchy proposed here can also be used to test whether 
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a set of Alice’s inputs is not jointly measurable. An interest¬ 
ing future direction would be to extend the present method to 
multipartite steering detection, in a quantum network scenario 
with some trusted and some untrusted parties | ^ ^. 
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Appendix: Supplemental material 

Proof of the equivalence between unsteerability and the 
existence of a separable model 

In the proof that follows we assume Bob’s Hilbert space to 
be arbitrary (continuous or discrete variable), while for sim¬ 
plicity we assume discrete outcomes for Alice. The general¬ 
ization of the proof to continuous outcomes will be immediate 
as we shall see. 

First, we recall that Bob’s assemblage is unsteerable 

by Alice’s inputs x = !,...,« (with corresponding outcomes 
Gx - 1 ,..., dx) iff it can be expressed as, 

q,\p(a\x, A)p^, Vx, a. (A.l) 

A 

The first part of the proof amounts to expressing (lA.ll) in a 
suitable form in terms of deterministic functions (i.e. the Kro- 
necker delta function) that will prove very helpful. The basic 
tool we utilize is the following identity, 

p(a\x, d) = z Sa,a, p(ax\x, d), (A.2) 

for a particular input x, while i5,j is the Kronecker delta. By 
inserting in (IA.2I) the identities, Za,for every 
input i + X, we get, 

pia\x, A) = z P(aill, A) ■ ■ ■ p(a„ln, A), (A.3) 

ai ...a,I 

where the summation over is implicitly included. Substi¬ 
tuting (IA.3I) back in the assemblage dA.ll) we get the desired 
expression, 

^a\x ~ ^ j (-^-4) 

a\...an 

where the unnormalized positive semidefinite operators 
> 0 correspond to, 

= ^q'i/7(ai|l,d)"-p(a„|n,d). (A.5) 

,1 

In the second part of the proof, we will show that one can al¬ 
ways define a separable model pab for Alice and Bob, and ap¬ 
propriate measurement operators for Alice, that can reproduce 
an arbitrary unsteerable assemblage (IA.41 i. Consider each in¬ 
put X of Alice, with outcomes Gx- 1,..., dx, to correspond to 
a fictitious observable (hermitian operator) Ax such that, 

^xWx)a — ^xWx^A- (.^-6) 

where the same outcomes Gx - I ,..., dx correspond to its real 
eigenvalues with |ai)A, ■.., |a„)Abeing the corresponding eigen¬ 
vectors. When Alice announces to Bob a pair (g, x), i.e. mea¬ 
sured input X and got outcome g, it will be considered equiv¬ 
alent as if she measured the observable Ax and got the eigen¬ 
value G as an outcome (with corresponding eigenvector |fl)). 


Note that such a correspondence x Ax can always be made, 
since the announced outcomes Gx always correspond to eigen¬ 
values of some observable. 

Next, Gssume that all the defined observables {Ai,...,A„) 
mutually commute, 

[A;„A,,-]=0, Vx9ix', (A.7) 


and, therefore, a joint basis exists that diagonalizes all Ax, Vx, 
simultaneously. We will show that these commuting observ¬ 
ables can reproduce the statistics of any unsteerable assem¬ 
blage by acting on a suitable separable state. Let us denote 
the vectors of this basis as {|ai ■ ■ • a„)), which sum to unity, 
2ai...a„ I'^i • ■ ■Gn)A{a\ ■ ■ - Gnl - 1, and are orthonormal, i.e., 

{g^...G^\gi...G,,} = ' ' ' da„,a',- (.^- 8 ) 

Due to the simultaneous diagonalization of every observable, 

it holds, Ax\gi ■ ■ ■ Gn)A — Gx\gi ■ ■ ■ Gn)A, Gx, X. 

If Pab is the shared state between Alice and Bob, when Al¬ 
ice measures input x Ax and announces output g, Bob’s 
(unnormalized) state conditioned on the pair (x, g), will be, 

^\X = T^rA[(Ma\x ® Ib)Pab], (A.9) 

where we defined the projectors onto the eigenstates of Ax 
with eigenvalue g, 

Ada\x ~ ^ ‘ ' ' ^n)A^^\ ' ' ' ^n\, (A. 10) 

satisfying, M^|^, = Ma\x and Zo Ma\x - 1, V x. Notice the sum¬ 
mation over the outcomes of the unannounced inputs, which is 
due to the inaccessibility of these degrees of freedom to Bob. 
Using the spectral decomposition of each Ax, we also get an 
expression for the observables, i.e., 

dx 

A, = (A. 11) 

a-l 


Now we will show the desired result that if pab is the fol¬ 
lowing separable state, 

Pab— 'y ] \^l, • ■ • , ^n}A{^l, • ■ • , ^n\® (A.12) 

Bob’s conditional state ( IA.9I ) will correspond to the unsteer¬ 
able assemblage (IA.4I) if Alice measures the commuting ob¬ 
servables defined in (lA.l 11 1 . We have. 


O'fp = ^U[{Ma\x ® 1b)PAb] 

(A. 13) 

“ ^ j da,ax^ai...a,i. 


matching exactly dAAl i. where we used the orthonormality of 
the states (IA. 8 I 1 and the property - dij of the Kronecker 
delta. 
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The generalization of the proof from discrete to continu¬ 
ous outcomes for Alice is straightforward, by replacing all 

CO 

summations with integrals, 2 ^ J da^ and the Kronecker 

—oo 

delta with the Dirac delta function, ^ 6{a - a^), which 
is a common practice when dealing with continuous Hilbert 
spaces. 


SDP, Dual and Optimal steering witnesses 


First, we will show that the problem (HJi can be expressed 
as an SDP 0371 Ihlll . and then derive its corresponding dual 
problem that will lead us to the optimal steering witnesses. 

Consider a square N xN (moment) matrix F with arbitrary 
elements r,^. Whether such a matrix is positive semidefinite, 
i.e. r > 0, is equivalent to whether its smallest eigenvalue, d*, 
is non-negative, i.e. d* > 0. The steering detection method 
outlined in the main Letter, boils down to finding the maxi¬ 
mized d* (name it, d™’‘) over all possible (complex, in gen¬ 
eral) values of the moment matrix’s elements {Fij} satisfying 
at the same time two types of constraints: 

(a) All the observable elements of F are constrained to be 
equal to the observable values from the steering test 1142 1. 

(b) Linear relations between the unobservable elements, im¬ 
posed by the commutativity constraint of Alice’s operators 
and the utilization of Bob’s operator algebra. 


The semidefinite program correspondin g to the problem de¬ 
scribed takes the following standard form 1371], 


imax 


max d 

-i.ILd 

subject to F - d 1 > 0 

Tr[FA,]=^', i = l,...,k 
Tr[FC,]=0, 7 = 1,...,/ 


(A. 14) 


known as the primal problem, the output of which will be 
dmax jjj (|a.141 i guarantees that the out¬ 

put of the SDP will be equal to the smallest eigenvalue of the 
given F. The second and third constraints correspond to the 
constraints (a) and (b) respectively, with suitably chosen ma¬ 
trices A, and Cj depending on the particular F, while k and / 
correspond to the total number of observable elements and lin¬ 
ear relations respectively. The values /?, are the ones obtained 
from the steering test, as explained in ll42ll . Concluding, steer¬ 
ing will be witnessed from the SDP (IA.14b if d™“ < 0. 


To obtain the dual of the SDP (IA.141 I. the solution of which 
will give us an upper bound on the quantity of interest d™‘*’‘. 


we start by writing the Lagrangian of this problem 116111 . 
i: = dH-Tr[Z-(F-dl)]-H 


+ ^ 71 * {bi - Tr[F • A,]) + 2 V* (o - Tr[F ■ Cj]) 

i=i j=i 

k r r /t / 

2 p*bi + d (1 - TrZ) + Tr F ■ Z - ^ p*Ai - ^ v)Cj 

7=1 y 

(A.15) 


(=1 


where the N x N hermitian matrix Z and the complex vari¬ 
ables {pi} and {vj} are the dual variables to the first, second 
and third (sets of) constraints in (IA.141 i respectively. If we 
consider the maximized value max £. over the primal vari- 

y,ir],7 

ables d, {F)y, it’s straightforward to see from (IA.15b that, 
maxX > d™“ H- Tr[Z ■ (F - d™“ 1)]. Therefore choosing 

Z > 0, and since F - d™™ 1 > 0 due to the first constraint in 
(lAT4l i. we find the following bound. 


maxi:>d”“, (A. 16) 

Our goal is to use max £. to get a good estimate for the figure 

of merit d™“, and in order for the bound (IA.16b not to be triv¬ 
ial (i.e. equal to infinity), £. should be bounded from above. 
We see that this occurs trivially if we set as constraints for the 
dual variables, 

TrZ=l (A. 17) 

k I 

+ (A. 18) 

/=l 7=1 

in addition to Z > 0. Imposing these constraints on X, the 
Lagrangian (lA.lSb optimized over the primal variables takes 
the simple form, 

k 

maxX = p*bi > d™“. (A. 19) 


Therefore we are lead to an alternative approach to bound the 
desired quantity d™“, by minimizing the left-hand side over 
the dual variables, for given {/?,), leading us to the following 
dual problem. 




- min 

iiiiUvjUZij] 

subject to 


1=1 

Z> 0 
TrZ = 1, 



(A.20) 


The output of the dual (IA.20b . y6*, is the tightest upper 
bound to the figure of merit d™“ ( IA.19b since optimal coef¬ 
ficients [pi] are found for the given observable values {bi}. A 
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negative value, y 6 * < 0 , is a sufficient condition for steerabil¬ 
ity, since it would imply that /i™“ < 0 , while a non-negative 
value > 0 is obtained for all unsteerable assemblages. 
Also, note that mere knowledge of the dual matrix Z (out¬ 
put of (IA.20b ) and the moment matrix F is enough to find 
since, Tr [FZ] = /?*, due to the second and third constraints in 
(IA.14b . To generalize this witness to any system, and therefore 
to arbitrary observations, consider arbitrary observable values 
{hi} + {bj} but keep the same coefficients {//,) as before. The 
following linear inequality, or steering witness, 

k 

Yjh:bi>0, (A.21) 

1=1 


detPz = 1 - <A2 ® Yf - (As ® Zf -RI>0, (A.24) 

detPs = l-<Ai®Z)2-<A3®Z)2-p2>0, (A.25) 

detP4 = 1 - <Ai ® - <A2 ® Yf - Pi > 0, (A.26) 

with, 

/ (Pi , P2, Ps) = (P3 <Ai ® Z> - P2 (A2 ® T) H- Pi (As ® Z)f 

-R\-Rl-Rl 

(A.27) 

where the matrix P, is obtained by F.^ by deleting its /-th row 
and column. Each of the conditions (IA.24b - (IA.26b leads to a 
steering criterion. For example. 


is satisfied by all unsteerable assemblages while a violation 
signals steering detection. For the particular {b,} the violation 
of (IA.21b is maximal since the coefficients (p,) are optimal 
for these particular values and non-optimal for any other, and 
therefore we refer to (IA.21b as the optimal steering witness for 
the values (I?,) = {bj}, obtained by particular measurements 
and assemblages. 

Finally, it is easy to verify that the primal problem is strictly 
feasible - i.e. there exists a F satisfying all the equality con¬ 
straints which is strictly positive definite. As such, strong du¬ 
ality holds for the primal and dual SDP problems, such that 
the optimal value of the primal 4™“ and the optimal value of 
the dual yS* are equal. 


Analytical derivation of non-linear steering criteria 


Consider the moment matrix F 51 ; (|5]l obtained by the set of 
measurements, >S = {ll®ll,Ai®A',A 2 ®T,As®Z), where the 
statistics of Alice’s unknown measurements Ai,A 2 ,As also 
originate from “spin”-measurements X, Y, Z. In the following 
derivation, only the algebra of Alice’s and Bob’s observables 
will matter, independently of their shared state pab- Applying 
the steps of the detection method, i.e. commutativity and the 
operator algebra on Bob’s side, the matrix (|5ll can be seen to 
get the simple form. 




1 (Ai®A) (A2®T) (As®Z>^ 
(Ai ® A) 1 iR\ iR 2 

(A 2 ®T) -/Pi 1 /Ps 

(As ® Z) —/P 2 —/Ps 1 > 


where the three free parameters P, are real, and equal to. 
Pi = (A 1 A 2 ® Z), P 2 = (AsAs ® X}, and Ps = —(AiAs ® T). 
Notice that the diagonal observable terms are equal to unity 
independently of the shared state, due to the fact that the Pauli 
operators, and the observables of Alice, take values +1, and 
therefore square to the identity. 

As explained in the main text, the necessary condition for 
unsteerability F^ > 0 implies the following conditions for its 
principal minors. 


detF^ = 1 -(Ai®A)2-(A2®T)^-(As®Z)2 
+ /(Pi,P2,Ps)>0, 


detP 2 > 0 1 - (A 2 ® Yf - (As ® Zf >rI>0, (A.28) 

and similarly for (IA.25I) . (IA.26I) . A violation of the last in¬ 
equality in (IA. 28 I 1 signals steering since there exist no assign¬ 
ment for the free parameters P, that can make (IA.28I) non¬ 
negative. When applied to the family of Werner states these 
criteria can be seen to detect steering for w > -^, which is 
a weaker detection than what the optimal witness (|6]) and the 
stronger non-linear criterion (|7]i can achieve. This is of course 
to be expected, since the former criteria only involve two mea¬ 
surement settings per site. 

The stronger non-linear criterion (|2l), based on three mea¬ 
surement settings, can be derived from (IA.231 I. where the con¬ 
tribution of the free parameters is grouped in the function 
/(Pi,P 2 ,Ps)- Our goal is to provide an upper bound for this 
function, say / < /max, and therefore limit its capability of 
making (IA.23I I positive for any given measurements. As a 
simple example of the logic behind, the analogous function in 
(IA.24b would be -P^ and is upper bounded by zero, as seen in 
the steering criterion ( IA.28I ). The maximum of /(Pi,P 2 ,Ps) 
can be seen to correspond to the following values for Pj, P 2 , 


duj - 


Therefore, 


Pf =Ps 
P* = -Ps 


(Ai ®A)(As ®Z) 
det P 2 + P 3 
(Ai ® A)(A2®y) 
det P 2 + R\ 


(A.29) 

(A.30) 


f{Rx,R2,R2)<f{R*,R*2,Rf) 

2 1 - (A, ® A)2 - (A 2 ® Yf - (As ® zf 


= -PS 


det P 2 + R\ 


(A.31) 


We employ this bound in (IA.23b and find that unsteerability 
of Bob’s assemblage implies. 


detF^i; > 0 ^ 

l-{Ai®Xf- (A2 ® Yf - (As ®Zf +f (Pf, P*, Ps) > 0 
« (1 - (Ai ® - (A2 ® Yf - (As ® Z) 2 ) , > 0 

(A.32) 


(A.23) 






































4 


Unsteerable assemblages necessarily satisfy detP 2 ^ 0 (see 
(IA.24b ). and therefore the last inequality of (IA.32b implies the 
desired non-linear criterion (|7]i, 

{Ai<S>Xf+ {A2<S>Yf+ {Ai«,Zf < 1. (A.33) 

Notice that for the expressions (IA.29b . (IA.30b we have as¬ 
sumed, KA 2 ® y)| < 1 and |(A 3 ®Z)| < 1. The cases where 
equality is attained in either (or both) inequalities should be 
treated separately, and it’s straightforward to see that in ev¬ 
ery single case the same condition (IA.33b is always obtained. 
Therefore, the validity of (IA.33b extends to the whole range 
of possible experimental outcomes. 

Optimal Witness for Lossy Single Photon state 

In this appendix we provide the optimal steering witness 
which certifies the steerability of the noisy single photon state. 
As described in the main text, we used the 11x11 moment 
matrix defined by the set of operators = {11®1 ,Ao®^b,Ao® 
Pb, Ai®^b, Ai®/7b, Aq® 1, Aj®l, l®qBPB, 1®Pb^b, 1® 


p|). First, note that moments of the form (A^. ® B) appearing 
in the moment matrix, with B an arbitrary string of length 2 
or more, are expected in general to be hard to measure experi¬ 
mentally. In the following we therefore assume these terms to 
be unobservable (and therefore treat them as free parameters 
in the moment matrix), and apply only the operator algebra of 
Bob to place linear relations between them. On the other hand, 
local moments of the form (1 ® Z?) can be measured efficiently 
by Bob, for example by estimating his local Wigner function 
or by using a linear optics scheme proposed by Shchukin and 
Vogel iSbll . and therefore we keep these moments as observ¬ 
able. The freedom that the method gives us to keep only 
those measurements that can be efficiently performed as ob¬ 
servable, highlights the flexibility of our approach to maintain 
experimental feasibility. Our ultimate goal is to provide an 
experimentally-friendly optimal steering witness. 

The code was implemented using cvx for matlab llb^l . with 
the optimal inequality extracted by solving the primal (lA. 141 i 
and dual (IA.201 i problems. The optimal inequality (IA.211 i for 
the noisy single photon state with q = 0.67 is given by 


P = 8.1657 - «Ao ® qB) + <Ai ® pb)) + 0.2508 «Ao ® ql) + <Ai ® pi)) - 0.3110 ((A^) + (A?)) 

-I- 0.3205 ((Aq ® ql) + <Aj ® pi)) + 0.3020 «Aq ® pi) + <Aj ® ql)) - 0.0001 «Aq ® qB) + <Aj ® pb)) 

+ 7.7217 «4) + {pi)) + 15.5451 {qlpl) - 31.0941 {{ql) + {pi)) - 2,1.0902,1 {qsps) > 0, (A.34) 


satisified by all unsteerable assemblages, with the state numer¬ 
ically achieving the violation p - -8.88 x 10 which is (in 
magnitude) far above the numerical precision. Smaller values 
of q still show a violation, with numerical evidence suggest¬ 
ing all ?7 > 2/3 demonstrate steering. The maximum violation 
of the inequality is /Zmax = -0.1556, achieved for q - \. 

Let us now comment on the experimental feasibility for 
the estimation of the witness (IA.34b . Most of the terms in 
Eq. (IA.34b can be efficiently measured by performing homo¬ 
dyne detection. The term that provides some extra difficulty 
in its measurement is the local fourth-order moment {qlpl) 
of Bob. As mentioned before, for the estimation of this term 
Bob could implement tomography on his local state, which 
doesn’t require conditioning on Alice’s outcomes. A more 
efficient approach that avoids tomography would be to use a 
scheme proposed by Shchukin and Vogel ifs^ . based on lin¬ 


r 


ear optics, that was designed to measure such local moments. 
A similar scheme was recently implemented by Avenhaus et 
al. 115711 . who managed to accurately measure moments of a 
single-mode up to eighth order. Therefore, we can safely con¬ 
clude that the proposed steering witness (IA.34b can be effi¬ 
ciently measured in the laboratory. 

Finally, let us note that the only terms which appear in 
the inequality are those which were considered observable in 
the moment matrix. However, observable terms of the form 
(A* ® B), which are experimentally demanding, were consid¬ 
ered unobservable, and as one would expect steering detec¬ 
tion weakens due to such relaxation. If on the other hand we 
consider all these experimentally demanding terms to be ob¬ 
servable, we find the same critical noise q > 2/3, with only 
the magnitude of the violation increasing (and the inequality 
containing the additional observable terms absent in (IA.34b '). 































